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Abstract. Building on work of Barker, Humpherys, Lafitte, Rudd, 
and Zumbrun in the shock wave case, we study stability of compressive, 
or "shock-like" , boundary layers of the isentropic compressible Navier- 
Stokes equations with 7-law pressure by a combination of asymptotic 
ODE estimates and numerical Evans function computations. Our re- 
sults indicate stability for 7 £ [1, 3] for all compressive boundary-layers, 
independent of amplitude, save for inflow layers in the characteristic 
limit (not treated). Expansive inflow boundary-layers have been shown 
to be stable for all amplitudes by Matsumura and Nishihara using en- 
ergy estimates. Besides the parameter of amplitude appearing in the 
shock case, the boundary-layer case features an additional parameter 
measuring displacement of the background profile, which greatly com- 
plicates the resulting case structure. Moreover, inflow boundary layers 
turn out to have quite delicate stability in both large-displacement and 
large-amplitude limits, necessitating the additional use of a mod-two 
stability index studied earlier by Serre and Zumbrun in order to decide 
stability. 
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1. Introduction 
Consider the isentropic compressible Navier-Stokes equations 

Pt + (pu)x = 0, 
(pu) t + (pu 2 ) x + p(p) x = u xx 

on the quarter-plane x, t > 0, where p > 0, u, p denote density, velocity, and 
pressure at spatial location x and time t, with 7-law pressure function 

(2) p(p) = a 0/ o 7 , a > 0, 7 > 1, 

and noncharacteristic constant "inflow" or "outflow" boundary conditions 

(3) (p,u)(0,t) = (p ,u ), u > 
or 

(4) n(0, t) =u u < 

as discussed in [251 EH EJ- The sign of the velocity at x = determines 
whether characteristics of the hyperbolic transport equation pt + up x = f 
enter the domain (considering / := pu x as a lower-order forcing term), and 
thus whether p(0, t) should be prescribed. The variable-coefficient parabolic 
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equation put — u xx = g requires prescription of u(0, t) in either case, with 
g ■= -p( u 2 /2) x -p(p) x - 

By comparison, the purely hyperbolic isentropic Euler equations 

^ Pt + (pu) x = 0, 

(pu) t + (pu 2 ) x +p(p) x = 

have characteristic speeds a = u ± \/p'(p), hence, depending on the values 
of (p, u)(0,t), may have one, two, or no characteristics entering the domain, 
hence require one, two, or no prescribed boundary values. In particular, 
there is a discrepancy between the number of prescribed boundary values 
for ([1]) and © in the case of mild inflow uq > small (two for ([I]), one 
for ©) or strong outflow uq < large (one for (TjQ), none for ©), indicat- 
ing the possibility of boundary layers, or asymptotically-constant stationary 
solutions of (P): 

(6) (p,u)(x,t) = (p,u)(x), lim (p,u)(z) = (p+,u + ). 

Indeed, existence of such solutions is straightforward to verify by direct com- 
putations on the (scalar) stationary-wave ODE; see [201 [23 13 [151 HH1 [9] or 
Section [2731 These may be either of "expansive" type, resembling rarefaction 
wave solutions on the whole line, or "compressive" type, resembling viscous 
shock solutions. 

A fundamental question is whether or not such boundary layer solutions 
are stable in the sense of PDE. For the expansive inflow case, it has been 
shown in [19] that all boundary layers are stable, independent of amplitude, 
by energy estimates similar to those used to prove the corresponding result 
for rarefactions on the whole line. Here, we concentrate on the complemen- 
tary, compressive case (though see discussion, Section [LTj) . 

Linearized and nonlinear stability of general (expansive or compressive) 
small- amplitude noncharacteristic boundary layers of (pQ) have been estab- 
lished in |191 [251 QUI HDJ- More generally, it has been shown in [101 [26] 
that linearized and nonlinear stability are equivalent to spectral stability, 
or nonexistence of nonstable (nonnegative real part) eigenvalues of the lin- 
earized operator about the layer, for boundary layers of arbitrary amplitude. 
However, up to now the spectral stability of large- amplitude compressive 
boundary layers has remained largely undetermined^ 

We resolve this question in the present paper, carrying out a systematic, 
global study classifying the stability of all possible compressive boundary- 
layer solutions of (pQ) . Our method of analysis is by a combination of asymp- 
totic ODE techniques and numerical Evans function computations, following 
a basic approach introduced recently in [HI [3j for the study of the closely re- 
lated shock wave case. Here, there are interesting complications associated 
with the richer class of boundary-layer solutions as compared to possible 



See, however, the investigations of [25] on stability index, or parity of the number of 
nonstable eigenvalues of the linearized operator about the layer. 
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shock solutions, the delicate stability properties of the inflow case, and, in 
the outflow case, the nonstandard eigenvalue problem arising from reduction 
to Lagrangian coordinates. 

Our conclusions are, for both inflow and outflow conditions, that com- 
pressive boundary layers that are uniformly noncharacteristic in a sense to 
be made precise later (specifically, u+ bounded away from 1, in the terminol- 
ogy of Section I2.3f) are unconditionally stable, independent of amplitude, on 
the range 7 £ [1,3] considered in our numerical computations. We show by 
energy estimates that outflow boundary layers are stable also in the charac- 
teristic limit. The omitted characteristic limit in the inflow case, analogous 
to the small-amplitude limit for the shock case should be treatable by the 
singular perturbation methods used in \22\ [7] to treat the small-amplitude 
shock case; however, we do not consider this case here. 

In the inflow case, our results, together with those of [19J, completely 
resolve the question of stability of isentropic (expansive or compressive) 
uniformly noncharacteristic boundary layers for 7 G [1,3], yielding uncon- 
ditional stability independent of amplitude or type. In the outflow case, we 
show stability of all compressive boundary layers without the assumption of 
uniform noncharacteristicity. 

1.1. Discussion and open problems. The small-amplitude results ob- 
tained in [I9j [TBI [23j [10] are of "general type", making little use of the 
specific structure of the equations. Essentially, they all require that the dif- 
ference between the boundary layer solution and its constant limit at |x| =00 
be small in L 1 ^ As pointed out in [10], this is the "gap lemma" regime in 
which standard asymptotic ODE estimates show that behavior is essentially 
governed by the limiting constant-coefficient equations at infinity, and thus 
stability may be concluded immediately from stability (computable by exact 
solution) of the constant layer identically equal to the limiting state. These 
methods do not suffice to treat either the (small-amplitude) characteristic 
limit or the large-amplitude case, which require more refined analyses. In 
particular, up to now, there was no analysis considering boundary layers 
approaching a full viscous shock profile, not even a profile of vanishingly 
small amplitude. Our analysis of this limit indicates why: the appearance 
of a small eigenvalue near zero prevents uniform estimates such as would be 
obtained by usual types of energy estimates. 

By contrast, the large- amplitude results obtained here and (for expansive 
layers) in [19] make use of the specific form of the equations. In particular, 
both analyses make use of the advantageous structure in Lagrangian coor- 
dinates. The possibility to work in Lagrangian coordinates was first pointed 
out by Matsumura-Nishihara [19] in the inflow case, for which the station- 
ary boundary transforms to a moving boundary with constant speed. Here 
we show how to convert the outflow problem also to Lagrangian coordinates, 

"^Alternatively, as in |19l 123] . the essentially equivalent condition that xv'(x) be small 
in L 1 . (For monotone profiles, J + °° \v — v+\dx — ± / +o °(?) — v+)dx — =p Jq +0 ° xv'dx.) 
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by converting the resulting variable-speed boundary problem to a constant- 
speed one with modified boundary condition. This trick seems of general 
use. In particular, it might be possible that the energy methods of [19] 
applied in this framework would yield unconditional stability of expansive 
boundary-layers, completing the analysis of the outflow case. Alternatively, 
this case could be attacked by the methods of the present paper. These are 
two further interesting direction for future investigation. 

In the outflow case, a further transformation to the "balanced flux form" 
introduced in |22j . in which the equations take the form of the integrated 
shock equations, allows us to establish stability in the characteristic limit 
by energy estimates like those of [18J in the shock case. The treatment of 
the characteristic inflow limit by the methods of \22\ [7] seems to be another 
extremely interesting direction for future study. 

Finally, we point to the extension of the present methods to full (non- 
isentropic) gas dynamics and multidimensions as the two outstanding open 
problems in this area. 

New features of the present analysis as compared to the shock case consid- 
ered in 0[12] are the presence of two parameters, strength and displacement, 
indexing possible boundary layers, vs. the single parameter of strength in 
the shock case, and the fact that the limiting equations in several asymp- 
totic regimes possess zero eigenvalues, making the limiting stability analysis 
much more delicate than in the shock case. The latter is seen, for example, 
in the limit as a compressive boundary layer approaches a full stationary 
shock solution, which we show to be spectrally equivalent to the situation of 
unintegrated shock equations on the whole line. As the equations on the line 
possess always a translational eigenvalue at A = 0, we may conclude exis- 
tence of a zero at A = for the limiting equations and thus a zero near A = 
as we approach this limit, which could be stable or unstable. Similarly, the 
Evans function in the inflow case is shown to converge in the large-strength 
limit to a function with a zero at A = 0, with the same conclusions; see 
Section [3] for further details. 

To deal with this latter circumstance, we find it necessary to make use 
also of topological information provided by the stability index of [2~Tj l8l [25]. 
a mod-two index counting the parity of the number of unstable eigenval- 
ues. Together with the information that there is at most one unstable zero, 
the parity information provided by the stability index is sufficient to de- 
termine whether an unstable zero does or does not occur. Remarkably, in 
the isentropic case we are able to compute explicitly the stability index for 
all parameter values, recovering results obtained by indirect argument in 
[25] , and thereby completing the stability analysis in the presence of a single 
possibly unstable zero. 
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2. Preliminaries 

We begin by carrying out a number of preliminary steps similar to those 
carried out in [3j [12] for the shock case, but complicated somewhat by the 
need to treat the boundary and its different conditions in the inflow and 
outflow case. 



2.1. Lagrangian formulation. The analyses of [12} [3] in the shock wave 
case were carried out in Lagrangian coordinates, which proved to be par- 
ticularly convenient. Our first step, therefore, is to convert the Eulerian 
formulation (TjQ) into Lagrangian coordinates similar to those of the shock 
case. However, standard Lagrangian coordinates in which the spatial vari- 
able x is constant on particle paths are not appropriate for the boundary- 
value problem with inflow/outflow. We therefore introduce instead "psuedo- 
Lagrangian" coordinates 

(7) x := / p(y,t)dy, t := t, 

Jo 

in which the physical boundary x = remains fixed at x = 0. 

Straightforward calculation reveals that in these coordinates ([I]) becomes 



vt - svx - Ux = a{t)v£ 

(8) 

U t - SUx +p(y)x - = <T[t)Ux 



on x > 0, where 

(9) s = - — , a(t) = m(t) - s, m(t) := -p(0,t)u(0,t) = -u(0,t)/v(0,t), 

so that m(t) is the negative of the momentum at the boundary x = x = 0. 
From now on, we drop the tilde, denoting x simply as x. 



2.1.1. Inflow case. For the inflow case, no > so we may prescribe two 
boundary conditions on (jHJ), namely 

(10) = v > 0, u\ x=0 = u >0 

where both uq,vq are constant. 



2.1.2. Outflow case. For the outflow case, uq < so we may prescribe only 
one boundary condition on ([8]), namely 

(11) u\ x=0 = u Q < 0. 

Thus v(0, t) is an unknown in the problem, which makes the analysis of the 
outflow case more subtle than that of the inflow case. 
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2.2. Rescaled coordinates. Our next step is to rescale the equations in 
such a way that coefficients remain bounded in the strong boundary-layer 
limit. Consider the change of variables 

(12) (x,t,v,u) — > (—esx,es 2 t,v/e,—u/(es)), 
where e is chosen so that 

(13) < v+ < v- = 1, 

where v+ is the limit as x — > +oo of the boundary layer (stationary solution) 
(v, u) under consideration and V- is the limit as x — > — oo of its continuation 
into x < as a solution of the standing-wave ODE (discussed in more detail 
just below). Under the rescaling (fl2j) . (j8]) becomes 

v t + v x - u x = a(t)v x , 
u t + u x + (av >) x = a(t)u x + — 

where a = ao£~ 7_1 s~ 2 , a = — u(0,t)/v(0,t) + 1, on respective domains 
x > (inflow case) a; < (outflow case) . 

2.3. Stationary boundary layers. Stationary boundary layers 

(v,u)(x,t) = (v,u)(x) 

of JIU) satisfy 

(a) v'-u' = 



(15) 



(6) n' + (av~^) = 

(c) (v,u)\ x =o = (v ,u ) 

(d) lim (v,u) = (v,u)±, 

x— >±oo 

where (d) is imposed at +oo in the inflow case, — oo in the outflow case and 
(imposing a = 0) uq = v . Using (fT5|) (a) we can reduce this to the study of 
the scalar ODE, 

(16) v' + (av^y=(j 

with the same boundary conditions at x = and x = ±oo as above. Taking 
the antiderivative of this equation yields 

(17) v' = H c (v) =v(v + av~~< + C), 

where C is a constant of integration. 

Noting that Tic 18 convex, we find that there are precisely two rest points 
of (|17h whenever boundary-layer profiles exist, except at the single parame- 
ter value on the boundary between existence and nonexistence of solutions, 
for which there is a degenerate rest point (double root of Tic)- Ignoring this 
degenerate case, we see that boundary layers terminating at rest point v + 
as x — > +oo must either continue backward into x < to terminate at a 
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second rest point V- as x — ► — oo, or else blow up to infinity as x — > — oo. 
The first case we shall call compressive, the second expansive. 

In the first case, the extended solution on the whole line may be recog- 
nized as a standing viscous shock wave; that is, for isentropic gas dynamics, 
compressive boundary layers are just restrictions to the half-line x > [resp. 
x < 0] of standing shock waves. In the second case, as discussed in [19], the 
boundary layers are somewhat analogous to rarefaction waves on the whole 
line. From here on, we concentrate exclusively on the compressive case. 

With the choice u_ = 1, we may carry out the integration of (I16p once 
more, this time as a definite integral from — oo to x, to obtain 

(18) v = H{v) = v{v -l + a{v-~< - 1)), 
where a is found by letting x — ► +oo, yielding 

no\ v+-l 7 l-v+ 

(19) a = -^I = ^T^I ; 

in particular, a ~ vT in the large boundary layer limit v+ — * 0. This is 
exactly the equation for viscous shock profiles considered in |12j . 

2.4. Eigenvalue equations. Linearizing (|14p about (£>,n), we obtain 

v(0,t) A , 



vt + v x - u x 



(m \ _(u E \ _v(0,t) 



(20) ^ + UX -\^ V )-\J)=— U 
(v,u)\x=o = (%(t),0) 

lim (v, u) = (0, 0) 

x— >+oo 

where = w(0), 

(21) h(v) = -v^ +1 + a( 7 - 1) + (a + 1)£ 7 
and {5, Ct denote perturbations of t), n. 

2.4.1. Inflow case. In the inflow case, u(0,t) = v(0,t) = 0, yielding 

At; + v x — u x = 

(22) fh(v) \ (u x 
An + Urr. — — 



An + u x — - , : v i , , 
on x > 0, with full Dirichlet conditions (v,u)\ x= q = (0,0). 

2.4.2. Outflow case. Letting U := (v,u) T , U := (v,u) T , and denoting by C 
the operator associated to the linearization about boundary-layer (y,u), 

(23) C := d x A(x) - d x B(x)d x , 
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(24) M*) = ( _ A(0 ) /0 , +1 "/ ) • BW = ( 2 °. 
we have t/j — £f7 = ^j^-U'(x), with associated eigenvalue equation 

(25) XU-CU= ?Ql*Lfr'(x), 

where U' = (v f , u'). 

To eliminate the nonstandard inhomogeneous term on the righthand side 
of ([25]), we introduce a "good unknown" (c.f. [3 © 031031) 

(26) U^u-X-^JM fr( x y 

Since CU' = by differentiation of the boundary-layer equation, the system 
expressed in the good unknown becomes simply 

(27) Ut-CU = in x < 0, 
or, equivalently, (|2"2"j) with boundary conditions 

C/U =0 = ^) (1 - A-V(O), -A~V(0)f 

(28) 

lim U = 0. 

Solving for u| x =o in terms of v| a;= o and recalling that f)' = v! by (|18p . we 
obtain finally 

-t)'(O) 

(29) u x=0 = a(AH x=0 , a(A) := — — . 

A — v (0) 

Remark 2.1. Problems ([25]) and (f27|) — ff22]) are evidently equivalent for all 
A ^ ; but are not equivalent for A = (for which the change of coordinates 
to good unknown becomes singular). For, U = U' by inspection is a solu- 
tion of (|27p . but is not a solution of (|25p . That is, we have introduced by 
this transformation a spurious eigenvalue at A = 0, which we shall have to 
account for later. 

2.5. Preliminary estimates. 

Proposition 2.2 ([3]). For each j > 1, < v + < 1/12 < v < 1, (|IBJ) 

has a unique (up to translation) monotone decreasing solution v decaying 
to endstates v± with a uniform exponential rate for v+ uniformly bounded 
away from u_ = 1. In particular, for < v+ < 1/12, 

(30a) \v(x) - v+\ < Ce'^ 1 x > 5, 

(30b) \v(x) -V-\< CeT^r 1 x<5 

where 5 is defined by v(5) = (v- + v+)/2. 
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Proof. Existence and monotonicity follow trivially by the fact that (|18f) is a 
scalar first-order ODE with convex righthand side. Exponential convergence 



as x — > +00 follows by H(v, v + ) = (v — v + ) [v— [ \_ v Ji j — ) ' wnence 

v — 7 < H jly^ < v — (1 — v + ) by 1 < < 7 for < x < 1. Exponential 
convergence as x — ► —00 follows by a similar, but more straightforward 
calculation, where, in the "centered" coordinate x := x — 5, the constants 
C > are uniform with respect to v + ,vo. See [3] for details. □ 

The following estimates are established in Appendices O and |Bj 



1- 



v ) 



Proposition 2.3. Nonstable eigenvalues X of (|22p . i.e., eigenvalues with 
nonnegative real part, are confined for any < w+ < 1 to the region 

(31) A:={A: Re(A) + |9fm(A)| < i^ + l) 2 }. 
for the inflow case, and to the region 

(32) A := {A : 5Re(A) + |9fm(A)| < max{^, 3 7 + §} 

2 8 

for the outflow case. 

2.6. Evans funct ion formulation. Setting w : — ~t~ v — u, we may 
express (f22j) as a first-order system 

(33) W' = A(x, X)W, 
where 



(34) A(x,\) 



A 


A ^ 











1. „.| 


M — V 


V 


m-xj 




V , / 



where 

(35) f(v) = v- v~^h(v) = 2v- 0(7 - l)v-"< - (a + 1), 
with /i as in (I21j) and a as in (jl9h . or, equivalently, 

(36) /(S) . M _ (7 _ 1)( 1^ )( ^_ ( 1^)„,_, 

Remark 2.4. T/ie coefficient matrix A may be recognized as a rescaled 
version of the coefficient matrix A appearing in the shock case [31 [12], with 

/l 0\ /l 
A= 1 \A 1 
\0 A/ \0 1/A, 

The choice of variables (w, u — v, v) T may be recognized as the modified flux 
form o/[22j, adapted to the hyperbolic-parabolic case. 
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Eigenvalues of (j22[) correspond to nontrivial solutions W for which the 
boundary conditions W(±oo) = are satisfied. Because A(x, A) as a func- 
tion of v is asymptotically constant in x, the behavior near x = ±00 of 
solutions of (|34l) is governed by the limiting constant-coefficient systems 

(37) W = A±(X)W, Afc(A) := A(±oo, A), 

from which we readily find on the (nonstable) domain 5RA > 0, A 7^ of in- 
terest that there is a one-dimensional unstable manifold W± (x) of solutions 
decaying at x = —00 and a two-dimensional stable manifold W^x) AW$~(x) 
of solutions decaying at x = +00, analytic in A, with asymptotic behavior 

(38) Wf{x, A) ~ e^ x)x V±{\) 

as x — ► ±00, where n±{\) and V^{\) are eigenvalues and associated ana- 
lytically chosen eigenvectors of the limiting coefficient matrices A±(X). A 
standard choice of eigenvectors Vp [EJ EJ EJ H3] , uniquely specifying Wj (up 
to constant factor) is obtained by Kato's ODE |15j . a linear, analytic ODE 
whose solution can be alternatively characterized by the property that there 
exist corresponding left eigenvectors Vp such that 

(39) (Vj ■ V ? ) ± = constant, (Vj ■ V J ) ± = 0, 
where " ' " denotes d/dX; for further discussion, see (T5j [HI [13]. 

2.6.1. Inflow case. In the inflow case, < x < +00, we define the Evans 
function D as the analytic function 

(40) An (A) := det«, W+, W+)\ x=0 , 

where Wj are as defined above, and W® is a solution satisfying the boundary 
conditions (v, u) = (0, 0) at x = 0, specifically, 

(41) W 1 °U=o = (1,0, Of. 

With this definition, eigenvalues of C correspond to zeroes of D both in 
location and multiplicity; moreover, the Evans function extends analytically 
to A = 0, i.e., to all of 3?A > 0. See [H EH [27] for further details. 
Equivalently, following |21[ [3] , we may express the Evans function as 

(42) D in (\) = (W°-W 1 + ) lx=Q , 

where W^~(x) spans the one-dimensional unstable manifold of solutions de- 
caying at x = +00 (necessarily orthogonal to the span of W^ix) and W%~(x)) 
of the adjoint eigenvalue ODE 

(43) W' = -A{x,X)*W. 

The simpler representation (|42l) is the one that we shall use here. 
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2.6.2. Outflow case. In the outflow case, — oo < x < 0, we define the Evans 
function as 

(44) Amt(A) := det(Wr,Wl W°) {x=0 , 

where is as defined above, and W® are a basis of solutions of ([33]) 
satisfying the boundary conditions ([29]) . specifically, 

(45) W 2 U =0 = (1,0, Of, W$\ x=0 = (0,- A _^ (Q) ,l^ T 
or, equivalently, as 

r>_^m = (wr . w9\ 

x=0' 



(46) D out (A) = (iVf 



where 



A 



(47) ^ = l°.- 1 .- A _ !i . ( o,. 

is the solution of the adjoint eigenvalue ODE dual to W® and W®. 

Remark 2.5. As discussed in Remark \2.1l D ont has a spurious zero at 
A = introduced by the coordinate change to "good unknown". 

3. Main results 
We can now state precisely our main results. 

3.1. The strong layer limit. Taking a formal limit as v + — > of the 

rescaled equations (|14p and recalling that a ~ vl, we obtain a limiting 
evolution equation 

v t + v x -u x = 0, 

(48) (u x \ 

u t + u x = — 

V V / x 

corresponding to a pressureless gas, or 7 = 0. 

The associated limiting profile equation v' = v(v — 1) has explicit solution 

n/ s l-tanh(^) 

(49) v°(x) = 2^^' 

v°(0) = 1 tanl ^( s / 2 ) = Vq -^ tfoQ limiting eigenvalue system is W' = A°(x, X)W, 

/OA A 

(50) A°(x,X) =00 A 

\v° v° f°(v°) — A, 

where f(v°) = 2v° - 1 = -tanh(^). 

Convergence of the profile and eigenvalue equations is uniform on any 
interval £>° > e > 0, or, equivalently, x — 5 < L, for L any positive constant, 
where the sequence of coefficient matrices is therefore a regular perturbation 
of its limit. Following [12], we call x < L+<5 the "regular region". For vq — > 
on the other hand, or x — > 00, the limit is less well-behaved, as may be seen 
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by the fact that df/dv ~ v^ 1 as v — > v + , a consequence of the appearance 
of y^-) in the expression (i36|) for /. Similarly, A(x, A) does not converge 
to Af (A) as i -> +00 with uniform exponential rate independent of 7, 
but rather as Cv~ 1 e~ x ^ 2 . As in the shock case, this makes problematic the 
treatment of x > L + 5. Following [T2] we call x > L + 5 the "singular 
region" . 

To put things in another way, the effects of pressure are not lost as v + — > 0, 
but rather pushed to x = +00, where they must be studied by a careful 
boundary-layer analysis. (Note: this is not a boundary-layer in the same 
sense as the background solution, nor is it a singular perturbation in the 
usual sense, at least as we have framed the problem here.) 

Remark 3.1. A significant difference from the shock case of [12] is the 

appearance of the second parameter vq that survives in the v+ — > limit. 

3.1.1. Inflow case. Observe that the limiting coefficient matrix 

/OA A \ 

(51) A° + (X) := A°(+oo,X) = A , 

\0 -l-A/ 

is nonhyperbolic (in ODE sense) for all A, having eigenvalues 0, 0, —1 — A; in 
particular, the stable manifold drops to dimension one in the limit v+ — > 0, 
and so the prescription of an associated Evans function is underdetermined. 

This difficulty is resolved by a careful boundary-layer analysis in [12J, 
determining a special "slow stable" mode 

y+±(i,o,o) T 

augmenting the "fast stable" mode 

y 3 :=(A/ / u)(A/^ + l),A/^,l) T 

associated with the single stable eigenvalue /i = —1 — A of A ,. This de- 
termines a limiting Evans function D® (X) by the prescription (|40p . (|38p of 
Section [2.61 or alternatively via (|42p as 

(52) Dl(\) = (W?°-W? + ) lx=0 , 

with W® + defined analogously as a solution of the adjoint limiting system 
lying asymptotically at x = +00 in direction 

(53) V x := (0,-1, A/Af 

orthogonal to the span of V2 and V3, where "~" denotes complex conju- 
gate, and Wf defined as the solution of the limiting eigenvalue equations 
satisfying boundary condition (14"T|) . i.e., (W° )i a;= o = (1,0, 0) T . 
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3.1.2. Outflow case. We have no such difficulties in the outflow case, since 
— oo) remains uniformly hyperbolic, and we may define a limiting 
Evans function D® ut directly by (pH]) , ([38]) , (j4"?| , at least so long as vq remains 
bounded from zero. (As perhaps already hinted by Remark 13.11 there are 
complications associated with the double limit (vo,v+) — > (0,0).) 

3.2. Analytical results. With the above definitions, we have the following 
main theorems characterizing the strong-layer limit i>+ — > as well as the 
limits Vq — > 0, 1. 

Theorem 3.2. For vq > n > and A in any compact subset of 3ftA > 0, 
-Di n (A) and D ou t(X) converge uniformly to D® n (\) and D® ut (\) as t>+ — ► 0. 

Theorem 3.3. For A in any compact subset of 5RA > and v + bounded 
from 1, D; n (A), appropriately renormalized by a nonvanishing analytic fac- 
tor, converges uniformly as v q —* 1 to the Evans function for the (uninte- 
grated) eigenvalue equations of the associated viscous shock wave connecting 
v_ = 1 to v + ; likewise, D® ut (\), appropriately renormalized, converges uni- 
formly as vo — > to the same limit for A uniformly bounded away from 
zero. 

By similar computations, we obtain also the following direct result. 

Theorem 3.4. Inflow boundary layers are stable for vq sufficiently small. 

We have also the following parity information, obtained by stability-index 
computations as in [25] 

Lemma 3.5 (Stability index). For any 7 > 1, vq, and v+, D[ n (0) 7^ ; hence 
the number of unstable roots of D in is even; on the other hand Df n (0) = 
and lim^oAaW = 0- Likewise, (D?J(0), D' out (0) ± 0, (D° out )'(0) ? 0, 
hence the number of nonzero unstable roots of Df n , D out , D® ut is even. 

Finally, we have the following auxiliary results established by energy es- 
timates in Appendices [Cl H2 [Ej and [FJ 

Proposition 3.6. The limiting Evans function D? is nonzero for A / on 
KeA > ; for all 1 > vq > 0. The limiting Evans function D® ut is nonzero 
for A ^ on $te\ > 0, for 1 > vq > v#, where ~ 0.0899 is determined by 
the functional equation v* = e -2 ^ 1- "*-* . 

Proposition 3.7. Compressive outflow boundary layers are stable for t> + 
sufficiently close to 1. 

Proposition 3.8 ([E]). Expansive inflow boundary layers are stable for all 
parameter values. 

Collecting information, we have the following analytical stability results. 



Indeed, these may be deduced from the results of [25] . taking account of the difference 
between Eulerian and Lagrangian coordinates. 
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Corollary 3.9. Forvo orv + sufficiently small, compressive inflow boundary 
layers are stable. For vo sufficiently small, v+ sufficiently close to 1, or 
vq > v* ~ .0899 and w+ sufficiently small, compressive outflow layers are 
stable. Expansive inflow boundary layers are stable for all parameter values. 

Stability of inflow boundary layers in the characteristic limit v+ — > 1 is 
not treated here, but should be treatable analytically by the asymptotic 
ODE methods used in (22J [7| to study the small-amplitude (characteristic) 
shock limit. This would be an interesting direction for future investigation. 
The characteristic limit is not accessible numerically, since the exponential 
decay rate of the background profile decays to zero as |1 — v+\, so that 
the numerical domain of integration needed to resolve the eigenvalue ODE 
becomes infinitely large as v+ — > 1. 

Remark 3.10. Stability in the noncharacteristic weak layer limit vq — > v+ 
[resp. 1] in the inflow [outflow] case, for w + bounded away from the strong 
and characteristic limits and 1 has already been established in [10\ 123] . 
Indeed, it is shown in [TO] that the Evans function converges to that for a 
constant solution, and this is a regular perturbation. 

Remark 3.11. Stability of Df n , D° ut may also be determined numerically, 
in particular in the region vq < i>* not covered by Proposition \3.6l 

3.3. Numerical results. The asymptotic results of Section [3.21 reduce the 
problem of (uniformly noncharacteristic, v+ bounded away from V- = 1) 
boundary layer stability to a bounded parameter range on which the Evans 
function may be efficiently computed numerically in a way that is uniformly 
well-conditioned; see [5]. Specifically, we may map a semicircle 

d{KA>0}n{|A| < 10} 

enclosing A for 7 G [1,3] by D® n , D® nt , D m , D out and compute the winding 
number of its image about the origin to determine the number of zeroes of 
the various Evans functions within the semicircle, and thus within A. For 
details of the numerical algorithm, see [2l[5]. 

In all cases, we obtain results consistent with stability; that is, a winding 
number of zero or one, depending on the situation. In the case of a single 
nonzero root, we know from our limiting analysis that this root may be quite 
near A = 0, making delicate the direct determination of its stability; how- 
ever, in this case we do not attempt to determine the stability numerically, 
but rely on the analytically computed stability index to conclude stability. 
See Section [6] for further details. 

3.4. Conclusions. As in the shock case [3| fT2] . our results indicate uncondi- 
tional stability of uniformly noncharacteristic boundary-layers for isentropic 
Navier-Stokes equations (and, for outflow layer, in the characteristic limit 
as well), despite the additional complexity of the boundary-layer case. How- 
ever, two additional comments are in order, perhaps related. First, we point 
out that the apparent symmetry of Theorem 13.31 in the vq — > outflow and 
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vq — ► 1 inflow limits is somewhat misleading. For, the limiting, shock Evans 
function possesses a single zero at A = 0, indicating that stability of inflow 
boundary layers is somewhat delicate as Vq — ► 1: specifically, they have an 
eigenvalue near zero, which, though stable, is (since vanishingly small in the 
shock limit) not "very" stable. Likewise, the limiting Evans function Df n as 
i>-t- — > possesses a zero at A = 0, with the same conclusions. 

By contrast, the Evans functions of outflow boundary layers possess a 
spurious zero at A = 0, so that convergence to the shock or strong-layer limit 
in this case implies the absence of any eigenvalues near zero, or "uniform" 
stability as v+ — * 0. In this sense, strong outflow boundary layers appear 
to be more stable than inflow boundary layers. One may make interesting 
comparisons to physical attempts to stabilize laminar flow along an air- or 
hydro-foil by suction (outflow) along the boundary. See, for example, the 
interesting treatise [24] . 

Second, we point out the result of instability obtained in [25] for inflow 
boundary-layers of the full (nonisentropic) ideal-gas equations for appropri- 
ate ratio of the coefficients of viscosity and heat conduction. This suggests 
that the small eigenvalues of the strong inflow-layer limit may in some cases 
perturb to the unstable side. It would be very interesting to make these 
connections more precise, as we hope to do in future work. 



4. Boundary-layer analysis 

Since the structure of (|34p is essentially the same as that of the shock 
case, we may follow exactly the treatment in [12] analyzing the flow of (|34p 
in the singular region x — ► +oo. As we shall need the details for further 
computations (specifically, the proof of Theorem l3.4p . we repeat the analysis 
here in full. 

Our starting point is the observation that 

/OA A \ 

(54) A(x, A) =00 A 

\v v f(v) - \) 

is approximately block upper-triangular for v sufficiently small, with diago- 
nal blocks ^ and (/(#) — A) that are uniformly spectrally separated 
on 3fteA > 0, as follows by 

(55) f(v) <v-l< -3/4. 

We exploit this structure by a judicious coordinate change converting (I34j) 
to a system in exact upper triangular form, for which the decoupled "slow" 
upper lefthand 2x2 block undergoes a regular perturbation that can be 
analyzed by standard tools introduced in [22]. Meanwhile, the fast, lower 
righthand lxl block, since scalar, may be solved exactly. 
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4.1. Preliminary transformation. We first block upper-triangularize by 
a static (constant) coordinate transformation the limiting matrix 

/OX X 

(56) A + = A(+oo, A) =00 A 

\v+ v+ /(«+)- A, 

at x = +oo using special block lower-triangular transformations 

< 57 > !)• L -=^' = (-M + ? 

where I denotes the 2x2 identity matrix and £ C lx2 is a 1 x 2 row 
vector. 

Lemma 4.1. On any compact subset of $te\ > 0, for each v + > swj^i- 
ciently small, there exists a unique 9+ = 6 + (v+, A) such that A + := L + A + R + 
is upper block-triangular, 



(58) i H 



A(J + u + 10+) Al 

/(«+) - A- Xv + 6 + l 



where J = ^jj anc? 11 = f ; satisfying a uniform bound 

(59) 1 0+| < C. 

Proof. Setting the 2 — 1 block of A + to zero, we obtain the matrix equation 

6+(aI - XJ) = -1 T + Xv + 9 + ie + , 
where a = f(v + ) — A, or, equivalently, the fixed-point equation 

(60) 6+ = (al- X J)' 1 (- l T + Xv + 6 + ld^ 



By det(al — XJ) = a 2 ^ 0, (al — AJ) _1 is uniformly bounded on compact 
subsets of 3fteA > (indeed, it is uniformly bounded on all of SJieA > 0), 
whence, for |A| bounded and v + sufficiently small, there exists a unique 
solution by the Contraction Mapping Theorem, which, moreover, satisfies 
(1591. □ 



4.2. Dynamic triangularization. Defining now Y := L + W and 

A(x, A) = L + A(x, X)R+(x, A) = 

x(j + v + ie + ) xi 

(v - v + )l T - v + (f(v) - f(v + ))9+ f(v) - A - Xv + Q + \ . 



we have converted (J34J) to an asymptotically block upper-triangular system 
(61) Y' = A(x, X)Y, 



L'R 
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with A + = A(+oo, A) as in (|58|) . Our next step is to choose a dynamic 
transformation of the same form 

(«) ^(|!). £=-r>-(^; 

converting (|6ip to an exactly block upper-triangular system, with G uni- 
formly exponentially decaying at x = +00: that is, a regular perturbation of 
the identity. 

Lemma 4.2. On any compact subset of JfteA > 0, for L sufficiently large 
and each v + > sufficiently small, there exists a unique G = G+(x,A, v + ) 
such that A := LA(x, X)R + L'R is upper block-triangular, 

rfi ox 7 A(j + u + i0+ + ie) ai \ 

[M) V /(£)-A-A0+l-AQlJ' 

and G(L) = 0, satisfying a uniform bound 

(64) |e(x,A,u+)| < Ce"^', 77 > 0, x > L, 
independent of the choice of L, v + . 

Proof. Setting the 2 — 1 block of A to zero and computing 

0\ (I 0\ _ ( o s 
-G' 0) {& I) ~ V-G' o„ 

we obtain the matrix equation 

(65) Q' - 8(aJ - A(J + u+10+)) = C + AGIO, 
where a(x) := f(v) — A — A-u + 6> + lL and the forcing term 

C := _(£ _ U+)1 T + _ /( u+ ))0 + 

by derivative estimate df/dv < CI) -1 together with the Mean Value Theorem 
is uniformly exponentially decaying: 

(66) ICI < C\v-v+\ < C 2 e^ x , 77 > 0. 

Initializing @(L) = 0, we obtain by Duhamel's Principle/Variation of 
Constants the representation (supressing the argument A) 

(67) &(x) = £ S^ x (C + X@ie)(y)dy, 

where S y ^ x is the solution operator for the homogeneous equation 
Q' - e(al - A(J + v+ld+)) = 0, 

or, explicitly, 

gV->x _ e fy a.(y)dy e -X(J+ v+ ie + )(x-y)^ 

For |A| bounded and v + sufficiently small, we have by matrix perturbation 
theory that the eigenvalues of — A( J + are small and the entries are 

bounded, hence 

i e -X(J+v + ie + )zi < q^z 
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for z > 0. Recalling the uniform spectral gap Sftea = /(«) — 3fteA < —1/2 for 
;ReA > 0, we thus have 

(68) \S y ^ x \ < Ce^'"^ 

for some C, r\ > 0. Combining ([66]) and ([68]) . we obtain 



(69) 



S y ^ x ((y)dy 



Defining O(x) =: 8(x)e (v/ 2 ) x anc [ recalling (|67p we thus have 

(70) 0(z) = / + e ( " /2)x ^ S y -* x e-^X9ie(y) dy, 

where / := e^/ 2 ^ J£ S y ^ x Q{y) dy is uniformly bounded, |/| < C3, and 
e (r]/2)x j% s y ^ x e~ m X6H9(y) dy is contractive with arbitrarily small contrac- 
tion constant e > in L°°[L, +00) for |#| < 2C3 for L sufficiently large, by 
the calculation 

( V /2)x f x s y ^ x e-™\0x^0i(y) -e^ l2)x f s y -" x e-^\e 2 ^e 2 {y) 



< 



( V /2)x / Ce - v (x-y) e - mdy | A |||^ 



y 2 ||oo max 1 1 t/j 1 1 00 

j 



< e ~(v/2)L 

= C 3 e- (T?/2)Z '|A| ||^i 



Ce -( v /2)(x-y) dy |A|||^ -^H^maxll^lU 

3 



72 00 max f,- oo- 

3 



It follows by the Contraction Mapping Principle that there exists a unique 
solution 9 of fixed point equation (fT0|) with |#(a;)| < 2C3 for x > L, or, 
equivalently (redefining the unspecified constant r/), (|64"|) . □ 

4.3. Fast /Slow dynamics. Making now the further change of coordinates 

z = Zy 

and computing 

(ZY)' = LY' + Z'y = (Za+ + L')Y, 

= (LA+R + L'R)Z, 

we find that we have converted (|6ip to a block-triangular system 

'\(J + V+19+ + 19) AIL \ 

/(£) - A - At>+6> + ll - Aeiy ' 

related to the original eigenvalue system (|34|) by 

I OA 
-e 1 r 



(71) Z' = AZ 



( 72 ) IY = LZ, i? := i? +J R = ( ^ J], I :— /?" 
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where 

(73) @ = e + v + e + . 

Since it is triangular, (|7ip may be solved completely if we can solve the 
component systems associated with its diagonal blocks. The fast system 

z' = (j(v) - A - A-M+l - XQljz 

associated to the lower righthand block features rapidly-varying coefficients. 
However, because it is scalar, it can be solved explicitly by exponentiation. 
The slow system 

(74) z' = \{J + v + 19 + + ie)z 

associated to the upper lefthand block, on the other hand, by (j64"|) . is an 
exponentially decaying perturbation of a constant-coefficient system 

(75) z = X(J + v + \9 + )z 

that can be explicitly solved by exponentiation, and thus can be well- 
estimated by comparison with (|75p . A rigorous version of this statement 
is given by the conjugation lemma of [20J: 

Proposition 4.3 ([20]). Let M(x,X) = M + (A) + Q(x,X), with M+ contin- 
uous in A and |G(x, A)| < Ce - ^, for A in some compact set A. Then, there 
exists a globally invertible matrix P(x, X) = I + Q(x, X) such that the coor- 
dinate change z = Pv converts the variable- coefficient ODE z' = M(x, X)z 
to a constant- coefficient equation 

v' = M + (X)v, 

satisfying for any L, < ry < r] a uniform bound 

(76) \Q{x,X)\ < C(L,f),r), max |(M + ) ij |, dimM + ) e -^' for x > L. 
Proof See [201E7], or Appendix C, [12] . □ 

By Proposition 14.3^ the solution operator for (|74p is given by 

(77) P(y, \) e HJ+v+w + (\,v + ))(x-y) p ( x ^ A )-i ; 

where P is a uniformly small perturbation of the identity for x > L and 
L > sufficiently large. 



5. Proof of the main theorems 
With these preparations, we turn now to the proofs of the main theorems. 
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5.1. Boundary estimate. We begin by recalling the following estimates 
established in [12] on W^{L + 5), that is, the value of the dual mode 
appearing in (|^2|) at the boundary x = L + 5 between regular and singular 
regions. For completeness, and because we shall need the details in further 
computations, we repeat the proof in full. 

Lemma 5.1 ([12]). For A on any compact subset of 3fteA > 0, and L > 
sufficiently large, with normalized as in [HJ [22j [3] , 

(78) \Wf(L + S)-V 1 \ < Ce-^ L 

as v + — * 0, uniformly in X, where C, rj > are independent of L and 

V X := (0,-l,A/ M ) T 
is the limiting direction vector (|53|) appearing in the definition of D® . 
Corollary 5.2 ([E]). Under the hypotheses of Lemma \5.1l 

(79) | W? + (L + 5) - Vi | < Ce~ nL 
and 

(80) \W+(L + 5)- W® + (L + S)\< Ce~ nL 

as — > 0, uniformly in A, where C , rj > are independent of L and W® + is 
the solution of the limiting adjoint eigenvalue system appearing in definition 
§2§ofD°. 



Proof of Lemma \5.1\ First, make the independent coordinate change x — ► 
x — 5 normalizing the background wave to match the shock-wave case. Mak- 
ing the dependent coordinate-change 

(81) Z := R*W, 

R as in (|72p . reduces the adjoint equation W = —A*W to block lower- 
triangular form, 

Z' = -A*Z = 

(82) f-\(j T + v ± ie + + iey o 
V — Ai T -f{v) + \ + \ v+ {e + i + Qiy 

with "~" denoting complex conjugate. 

Denoting by a suitably normalized element of the one-dimensional 
(slow) stable subspace of —A*, we find readily (see [12J for further discussion) 
that, without loss of generality, 

(83) V+^(0,l,H7 + X)- l f 

as v+ — > 0, while the associated eigenvalue flf — ► 0, uniformly for A on 
an compact subset of 5?eA > 0. The dual mode Z^ = R*W^ is uniquely 
determined by the property that it is asymptotic as x — ► +oo to the corre- 
sponding constant-coefficient solution e^ 1 V : (the standard normalization 

of PEE! E]). 
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By lower block-triangular form (|82|) . the equations for the slow variable 
z T := (Z\, Z2) decouples as a slow system 



(84) z' = -[X{J + v+\6 + + 16)J z 
dual to (|74p . with solution operator 

(85) P*(x, \)- l e -~ x{ - J+v ^ e +^-y)p{y, \y 

dual to ([77]). i.e. (fixing y = L, say), solutions of general form 

(86) z{\, x) = P*(x, A)- 1 e -^( J +^+ M +)*)(^) ? }, 

v G C 2 arbitrary. 
Denoting by 

Z+(L) := irW+(L), 

therefore, the unique (up to constant factor) decaying solution at +00, and 
vf := ((V±~)i, (V r 1 + )2) T > w e thus have evidently 

5+(x, A) = P*(x, X)' 1 e~^ J+v+19+ ^^ x vf , 
which, as v + — * 0, is uniformly bounded by 

(87) \zf(x,X)\<Ce ex 



for arbitrarily small e > and, by (|83p . converges for x less than or equal 
to X — 5 for any fixed X simply to 

lim zt(x,\) = P*(x,X)-\0,l) T . 

Defining by q := (Zf)% the fast coordinate of Z^~, we have, by (|82p . 
<? + (f(v) - A - (A?; + 6> + lL + A01)*)g = Al T z+, 



whence, by Duhamel's principle, any decaying solution is given by 
q(x,X) = [ +C ° e f v a{z ^ v+)dz Xl T z+(y)dy, 

J X 

where 

a(y, A, v+) := - (/(*>) - A — (Xv + + l + A01)*) . 

Recalling, for 3fteA > 0, that 5ftea > 1/2, combining (|87p and (|88p . and noting 
that a converges uniformly on y < Y as t> + — > for any Y > to 

a Q (y,A) := -/ («) + A + (A9 l)* 
= (1 + A) + 0(e-^) 
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we obtain by the Lebesgue Dominated Convergence Theorem that 

Jy ao{z ' X)dz A1 T (0, if dy 

e (l+X)(L-y)+J y L 0(e-*')dz dy 

= A(l + A)- 1 (l + 0(e- r ' L )). 
Recalling, finally, (I88j) . and the fact that 

\P-Id\(L,X), \R-Id\(L,X) < Ce~ nL 
for v+ sufficiently small, we obtain (|78p as claimed. □ 

Proof of Corollary \5.SX Again, make the coordinate change x — > x — 5 nor- 
malizing the background wave to match the shock-wave case. Applying 
Proposition 14.31 to the limiting adjoint system 

/ 0_ \ 
W' = —(A°)*W = j —A \W + 0(e-' nx )W, 
\-l -1 1 + Xj 

we find that, up to an Id + 0(e~ vx ) coordinate change, is given by 

the exact solution W = V\ of the limiting, constant-coefficient system 

/ 0_ \ 
W' = -{A°)*W = -A \W. 

\-l -1 1 + A/ 

This yields immediately (f79|l . which, together with (f78l) . yields (l80]l . □ 

5.2. Convergence to D°. The rest of our analysis is standard. 

Lemma 5.3. On x < L — 5 for any fixed L > 0, there exists a coordinate- 
change W = TZ conjugating (I34D to the limiting equations (I50p . T = 
T(x, A,v+), satisfying a uniform bound 

(89) \T - Id\ < C(L)v+ 

for all v + > sufficiently small. 

Proof. Make the coordinate change x — > x — 5 normalizing the background 
profile. For x £ (— oo,0], this is a consequence of the Convergence Lemma 
of [22], a variation on Proposition 14.31 together with uniform convergence 
of the profile and eigenvalue equations. For x £ [0,L], it is essentially 
continuous dependence; more precisely, observing that \A — A°\ < Ci(L)v + 
for x G [0, L], setting S := T — Id, and writing the homological equation 
expressing conjugacy of (|34j) and (|50l) . we obtain 

S'-(AS-SA°) = (A — A ), 
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which, considered as an inhomogeneous linear matrix-valued equation, yields 
an exponential growth bound 

S(x) < e Cx {S{0) + CT 1 Cx(L)v + ) 

for some C > 0, giving the result. □ 

Proof of Theorem \3.£% inflow case. Make the coordinate change x — > x — 5 
normalizing the background profile. Lemma 15.31 together with convergence 
as v+ —* of the unstable subspace of A- to the unstable subspace of A®_ 
at the same rate 0(v+) (as follows by spectral separation of the unstable 
eigenvalue of A and standard matrix perturbation theory) yields 

(90) \W?(0, A) - W?°(0, A) | < C(L)v+. 

Likewise, Lemma 15.31 gives 

\W+(0,\)-W? + (0,X)\ <C(L)v+\W+(0,X)\ 

+ \S^°\\W+(L,X)-W° + (L,X)\, 

where Sq~~* x denotes the solution operator of the limiting adjoint eigenvalue 
equation W' = — (A°)*W. Applying Proposition 14.31 to the limiting system, 
we obtain 

|<S L ^°| < C 2 e~ A °+ L < C 2 L\X\ 

by direct computation of e~ A+L , where C 2 is independent of L > 0. Together 
with dSDJ) and this gives 

\W+(0,X)-W° + (0,X)\ < C(L)v + \W+(0,X)\ + L\X\C 2 Ce^ L , 

hence, for |A| bounded, 

|^! + (0,A) - W?+(0,X)\ < C 3 (L)v + \W° + (0,X)\+LC A e-i L 

< C 5 (L)v + + LC^ L . 

Taking first L — ► oo and then v+ — ► 0, we obtain therefore convergence of 
Wf(0, A) and W^{0, A) to W? + (0, A) and ^+(0, A), yielding the result by 
definitions (gSJ) and ([52]). □ 

Proof of Theorem \3.SX outflow case. Straightforward, following the previous 
argument in the regular region only. □ 

5.3. Convergence to the shock case. 



Proof of Theorem \3.4\ inflow case. First make the coordinate change x — ► 
x — 5 normalizing the background profile location to that of the shock wave 
case, where 5 — ► +oo as Vq — ► 1. By standard duality properties, 

An = W? ■ W+\ x=Xo 

is independent of xq, so we may evaluate at x = as in the shock case. 
Denote by W-f, the corresponding modes in the shock case, and 

T> = WT-W?\ X=0 
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the resulting Evans function. 

Noting that VV+ and W\ are asymptotic to the unique stable mode at 
+00 of the (same) adjoint eigenvalue equation, but with translated decay- 
rates, we see immediately that VV^ = W+e - ^^ 1 . W® on the other hand is 
initialized at x = — 5 (in the new coordinates x = x — 5) as 

W 1 °(-5) = (1,0,0) T , 

whereas W{~ is the unique unstable mode at —00 decaying as e^ 1 X V{~ , where 
is the unstable right eigenvector of 

/OA A \ 
A- = A 

V 1 /a) -V 

Denote by the associated dual unstable left eigenvector and 

n r -.= vr(Vf) T 

the eigenprojection onto the stable vector . By direct computation, 
V,- =c(X)(l,l + X/^,^) T , c(A)/0, 

yielding 

(93) HyWf =: /3(A) = c(A) / 

for 5RA > 0, on which UfiJ > 0. 

Once we know (|93p . we may finish by a standard argument, concluding by 
exponential attraction in the positive x-direction of the unstable mode that 
other modes decay exponentially as x — > 0, leaving the contribution from 
/3(A) Vj~ plus a negligible 0(e - ' 7 ' 5 ) error, ?? > 0, from which we may conclude 

that W-fl^o ~ /3 -1 e~ 5At i WPl^o. Collecting information, we find that 

25(A) = /?(A)- 1 e- <5( ^+^ )(A) An(A) +0(e-" 5 ), 
r] > 0, yielding the claimed convergence as Vq — ► 1, 5 — > +00. □ 

Proof of Theorem \3.4\ outflow case. For A uniformly bounded from zero, 
W® = (0, —1, — A/(A — v'(0))) T converges uniformly as vq —* to 

(0,-1, -i) T , 

whereas the shock Evans function 25 is initiated by VV^ proportional to 

V+ = (0,-l,-l-A) T 

agreeing in the first two coordinates with W®. By the boundary-layer anal- 
ysis of Section [57T] the backward (i.e., decreasing x) evolution of the adjoint 
eigenvalue ODE reduces in the asymptotic limit v + — > (forced by vq — > 0) 
to a decoupled slow flow 

w = I Jw, w eC 
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in the first two coordinates, driving an exponentially slaved fast flow in the 
third coordinate. From this, we may conclude that solutions agreeing in the 
first two coordinates converge exponentially as x decreases. Performing an 
appropriate normalization, as in the inflow case just treated, we thus obtain 
the result. We omit the details, which follow what has already been done in 
previous cases. □ 

5.4. The stability index. Following [251 [10], we note that D m (\) is real 
for real A, and nonvanishing for real A sufficiently large, hence sgnD; n (+oo) 
is well-defined and constant on the entire (connected) parameter range. The 
number of roots of D m on 3? A > is therefore even or odd depending on the 
stability index 

sgn[D in (0)D in (+oo)]. 

Similarly, recalling that D out (0) = 0, we find that the number of roots of 
D out on 3? A > is even or odd depending on 

sgn[D^ ut (0)D out (+cx))]. 
Proof of Lemma \3. 5[ inflow case. Examining the adjoint equation at A = 0, 

W' = -A*W, -A*(x,0) = 





—f(v+) > 0, we find by explicit computation that the only solutions that are 
bounded as x — > +oo are the constant solutions W = (a,b,0) T . Taking the 
limit V^ + (0) as A — > + along the real axis of the unique stable eigenvector 
of —A+(\), we find (see, e.g., [27]) that it lies in the direction (1, 2 + 0) T , 
where > is the positive characteristic speed of the hyperbolic convection 

matrix (_ h ( v ^y v i+i > V i = c ( u o, v + )(l, 2 + a+, 0) T , c(v ,v + ) / 

0. Thus, .Din(O) = V{~ • (1,0, 0) r = c(vq,v + ) ^ as claimed. On the other 
hand, the same computation carried out for Df n (0) yields -0^(0) = 0. (Note: 

— 1/2 

a,j ~ v + — > +cxd as v+ — > 0.) Similarly, as vq — > 0, 

A° n (A)^(l,0,0) T -(0,l,*) T = 0. 

Finally, note -Di n (0) 7^ implies that the stability index, since continuously 
varying so long as it doesn't vanish and taking discrete values ±1, must be 
constant on the connected set of parameter values. Since inflow boundary 
layers are known to be stable on some part of the parameter regime by 
energy estimates (Theorem 13. 4p . we may conclude that the stability index 
is identically one and therefore there are an even number of unstable roots 
for all 1 > v > v + > 0. 

To establish that (-D?J'(0) / 0, we compute 

Df n '(0) = (dxW?°) ■ wl + + w?° ■ (d x w^+). 
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Since W® = (1,0,0) is independent of A, we need only show that the first 
component of d\ W® + is nonzero. Note that d\Wi + solves the limiting 
adjoint variational equations 

(d A W? + )'(0) + (A°)*(x,0)d x W° + = b(x) 

with 

/0 v° \ ( 

30° - 1 

By (f53j) . and the fact that d x p>i + = 0, 9 A Wi (a;) is chosen so that asymp- 
totically at x = +00 it lies in the direction of d\V\ = (0,0,-1). Set 
d x W? + = (d x W^ 1 ,d x W^ 2 ,d x W^ 3 ) T . Then the third component solves 

{d x W°+)' + v°d x W?+ = 63 := 3f}° - ^ - 1 
which has solution 








v° \ 




G 





v° , 


b(x) = 1 







/V)/ 





d x W?+(x) = d A <+(+ooMx) - <p(x) / <p-\y)h(v)dy 

J X 

where 

^(x) = e fr 

Integrating the equation for the first component of d x W® + yields 

rOC 

d x W?+(x) = d A <+(+oo) + / d x W°+(y)dy 

J X 

roo 

= d x W?+(+oo) + d x W°+(+oo) / v°(y)tp(y)dy 

J X 

poo / poo \ 

-J yf(y) J ^ 1 (z)b 3 (z)dzj dy. 

Using the condition d x W° + (+oo) = (0,0, -1) T we have 5 A W 1 04 L (+oo) = 
0,d x W°+(+oo) = -1 so that 

f'OQ f'OQ / POO \ 

d x W^\ x=0 = - J o v°(y)ip(y)dy- J ^(y) J <p-\z)b 3 dz) dy. 

Finally, note that by using (|49p we have 63 = 1 — tanh(^-^) so that for all 
x > 0, tp{x),bz{x) > which implies 

A°„ , (o) = 3 A ^ 1 0+ 1 |, =0 /o. 

□ 
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Remark 5.4. The result L>; n (0) 7^ at first sight appears to contradict that 
of Theorem \3.3\ since T>(0) = for the shock wave case. This apparent 
contradiction is explained by the fact that the normalizing factor e~ s ^ 1 + ^ > 
is exponentially decaying in 5 for A = ; since /i^"(0) = 0, while > 0. 
Recalling that 5 — * +00 as vq — * 1, we recover the result of Theorem [X 



Proof of Lemma \3.5[ outflow case. Similarly, we compute 

D f wst {0) = d\Wr-W? 1 

where d\Wi\\=o satisfies the variational equation Ld\U±(0) = U± = U', 
or, written as a first-order system, 

/ u x \ /0 

(0AWf)' -A(x,0)d\W{- = v x , A(x,0)= 

\-VxJ \v v f(v) 

which may be solved exactly for the unique solution decaying at — 00 of 
Wf(0) = I I , (0AWf)(O) = f v - v- 

Recalling from gZD that W?(A) = (0, -1, -A/(A - w'(0))) T , hence 

W?(0) = (0,-1, 0) T , 5 A ^°(0) = (0,0,1/^'(0)) T , 
we thus find that 

£>out(0) = d x Wr(0) ■ W?(0) + Wf(0) • W°(0) 
= -(t)(0)-l) + l = 2-uo^O 

as claimed. The proof that (D® nt )'(0i) ^ goes similarly. 

Finally, as in the proof of the inflow case, we note that nonvanishing 
implies that the stability index is constant across the entire (connected) pa- 
rameter range, hence we may conclude that it is identically one by existence 
of a stable case (Corollary 13.91) . and therefore that the number of nonzero 
unstable roots is even, as claimed. □ 

5.5. Stability in the shock limit. 

Proof of Corollarv \3.9c inflow case. By Proposition 13. 61 we find that D m has 
at most a single zero in IRA > 0. However, by our stability index results, 
Theorem 13.51 the number of eigenvalues in 5RA > is even. Thus, it must 
be zero, giving the result. □ 

Proof of Corollary \3.tk outflow case. By Theorem l3.3( D oxlt , suitably renor- 
malized, converges as v — > to the Evans function for the (unintegrated) 
shock wave case. But, the shock Evans function by the results of O [12] 
has just a single zero at A = on 3? A > 0, already accounted for in D ou t 
by the spurious root at A = introduced by recoordinatization to "good 
unknown" . □ 
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5.6. Stability for small vq. Finally, we treat the remaining, "corner case" 
as v+, vo simultaneously approach zero. The fact (Lemma 13. 5ft that 

lim lim An (A) = 

shows that this limit is quite delicate; indeed, this is the most delicate part 
of our analysis. 



Proof of Theorem \3.4\ inflow case. Consider again the adjoint system 

/0 v 

W' = -A*(x,X)W, A*(x,\)= ( A v 
By the boundary analysis of Section 15.11 

v -/(«)+ A / 

where a := , and fi is the unique stable eigenvalue of satisfying 
(by matrix perturbation calculation) 

~^=~X{vT + 0{v + )) 

1/2 

and thus a = v + + 0(v+) as vq — > (hence v + — > 0) on bounded subsets 
of 5RA > 0. Combining these expansions, we have 

W x {+oo) = v'fil + o(l)), W 3 = ~ A (1 + (1)) 

-/(«) + A 

for vq sufficiently small. 

From the W\ equation W = vW%, we thus obtain 

r+oo 

Wi(0) = Wi(+oo) - / vW 3 (y)dy 
Jo 

Observing, finally, that, for 5RA > 0, the ratio of real to imaginary parts of 
-f($)+\ (v) ^ s uniformly positive, we find that 3f?VVi(0) ^ for sufficiently 
small, which yields nonvanishing of -D;n(A) on 3?A > as claimed. □ 

6. Numerical computations 

In this section, we show, through a systematic numerical Evans function 
study, that there are no unstable eigenvalues for 

(7,«+) G [1,3] x (0,1], 

in either inflow or outflow cases. As defined in Section 12.61 the Evans func- 
tion is analytic in the right-half plane and reports a value of zero precisely 
at the eigenvalues of the linearized operator ([20]) . Hence we can use the 
argument principle to determine if there are any unstable eigenvalues for 
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(c) (d) 

Figure 1. Typical examples of the inflow case, showing con- 
vergence to the limiting Evans function as v+ — > for a 
monatomic gas, 7 = 5/3, with (a) vq = 0.1, (b) vo = 0.2, 
(c) vo = 0.4, and (d) vo = 0.7. The contours depicted, going 
from inner to outer, are images of the semicircle 4> under D 
for v + = le— 2, le— 3, le— 4, le— 5, le— 6, with the outer-most 
contour given by the image of cj> under D°, that is, when 
v + = 0. Each contour consists of 60 points in A. 

this system. Our approach closely follows that of j3j [12] for the shock case 
with only two major differences. First, our shooting algorithm is only one 
sided as we have the boundary conditions (|4ip and (I47|) for the inflow and 
outflow cases, respectfully. Second, we "correct" for the displacement in the 
boundary layer when vq ~ 1 in the inflow case and Vq ~ in the outflow 
case so that the Evans function converges to the shock case as studied in 
[3"1 [T2] (see discussion in Section HOj) . 

The profiles were generated using Matlab's bvp4c routine, which is an 
adaptive Lobatto quadrature scheme. The shooting portion of the Evans 
function computation was performed using Matlab's ode45 package, which 
is the standard 4th order adaptive Runge-Kutta-Fehlberg method (RKF45). 
The error tolerances for both the profiles and the shooting were set to 
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(c) (d) 

Figure 2. Typical examples of the outflow case, showing 
convergence to the limiting Evans function as v + — » for a 
monatomic gas, 7 = 5/3, with (a) vo = 0.2, (b) vq = 0.4, (c) 
vo = 0.6, and (d) vq = 0.8. The contours depicted are images 
of the semicircle <j> under D for v+ = le— 2, le— 3, le— 4, le— 
5, le — 6, and the limiting case v+ = 0. Interestingly the 
contours are essentially (visually) indistinguishable in this 
parameter range. Each contour consists of 60 points in A 

AbsTol=le-6 and RelTol=le-8. We remark that Kato's ODE (see Sec- 
tion 12.61 and [151 E] for details) is used to analytically choose the initial 
eigenbasis for the stable/unstable manifolds at the numerical values of infin- 
ity at L = ±18. Finally in SectionEH we carry out a numerical convergence 
study similar to that in [3]. 

6.1. Winding number computations. The high-frequency estimates in 
Proposition 12.31 restrict the set of admissible unstable eigenvalues to a fixed 
compact triangle A in the right-half plane (see (I3ip and (1321) for the inflow 
and outflow cases, respectively). We reiterate the remarkable property that 
A does not depend on the choice of u+ or vq. Hence, to demonstrate stability 
for a given 7, v+ and vq, it suffices to show that the winding number of the 
Evans function along a contour containing A is zero. Note that in our region 
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(a) (b) 

Figure 3. Typical examples of the Evans function evaluated 
along the positive real axis. The (a) inflow case is computed 
for vq = 0.7 and vq = and (b) the outflow case is computed 
for vq = 0.3 and v+ = 0.001. Not the transversality at the 
origin in both cases. Both graphs consist of 50 points in A. 

of interest, 7 6 [1,3], the semi-circular contour given by 

<f>:=d{{\ I 5Re> 0} U {A | |A| < 10}), 

contains A in both the inflow and outflow cases. Hence, for consistency we 
use this same semicircle for all of our winding number computations. 

A remarkable feature of the Evans function for this system, and one that is 
shared with the shock case in [31 [12] , is that the Evans function has limiting 
behavior as the amplitude increases, Section [3T2l For the inflow case, we see 
in Figure [H the mapping of the contour <fi for the monatomic case (7 = 5/3), 
for several different choices of Vq, as v+ — ► 0. We remark that the winding 
numbers for < v + < 1 are all zero, and the limiting contour touches zero 
due to the emergence of a zero root in the limit. Note that the limiting 
case contains the contours of all other amplitudes. Hence, we have spectral 
stability for all amplitudes. 

The outflow case likewise has a limiting behavior, however, all contours 
cross through zero due to the eigenvalue at the origin. Nonetheless, since 
the contours only wind around once, we can likewise conclude that these 
profiles are spectrally stable. We remark that the outflow case converges 
to the limiting case faster than the inflow case as is clear from Figure [2j 
Indeed, v+ = le—2 and the limiting case v+ = 0, as well as all of the values 
of i>+ in between, are virtually indistinguishable. 

In our study, we systematically varied vq in the interval [.01, .99] and 
took the v + — > limit at each step, starting from a v+ = .9 (or some 
other appropriate value, for example when vq < .9) on the small-amplitude 
end and decreased v + steadily to 10~ fc for k = 1,2,3, ... ,6, followed by 
evaluation at v + = 0. For both inflow and outflow cases, over 2000 contours 
were computed. We remark that in the v + — > limit, the system becomes 



STABILITY OF BOUNDARY LAYERS 



33 




(a) (b) 

Figure 4. Shock limit for (a) inflow and (b) outflow cases, 
both for 7 = 5/3. Note that the images look very similar to 
those of PES]. 



pressureless, and thus all of the contours in the large-amplitude limit look 
the same regardless of the value of 7 chosen. 

6.2. Nonexistence of unstable real eigenvalues. As an additional ver- 
ification of stability, we computed the Evans function along the unstable 
real axis on the interval [0, 15] for varying parameters to show that there 
are no real unstable eigenvalues. Since the Evans function has a root at 
the origin in the limiting system for the inflow case, and for all values of 
v + in the outflow case, we can perform in these cases a sort of numerical 
stability index analysis to verify that the Evans function cuts transversely 
through the origin and is otherwise nonzero, indicating that there are no 
unstable real eigenvalues as expected. In Figure EJ we see a typical example 
of (a) the inflow and (b) outflow cases. Note that in both images, the Evans 
function cuts transversally through the origin and is otherwise nonzero as A 
increases. 

6.3. The shock limit. When vq is far from the midpoint (1 — v+)/2 of 
the end states, the the Evans function of the boundary layer is similar to 
the Evans function of the shock case evaluated at the displacement point 
xq. Hence, when we compute the boundary layer Evans function near the 
shock limits, vq ~ 1 for the inflow case and v ~ for the outflow case, we 
multiply for the correction factor c(A) so that our output looks close to that 
of the shock case studied in [31 [12] . The correction factors are 

C (A) = e (-M + -A-)*o 

for the inflow case and 

C (A) = e (-A + -/0*o ; 

for the outflow case, where fi~ is the growth mode of A__(A) and n + is the 
decay mode of A + (X). In Figure HJ we see that these highly displaced profiles 
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Inflow Case 


L 


7 = 1.2 


7 = 1.4 


7 = 1.666 


7 = 2.0 


7 = 2.5 


7 = 3.0 


8 


7.8(-l) 


8.4(-l) 


9.2(-l) 


1.0(0) 


1.2(0) 


1.3(0) 


10 


1.4(-1) 


1.2(-1) 


9.2(-2) 


6.8(-2) 


4.4(-2) 


2.8(-2) 


12 


1.4(-2) 


7.9(-3) 


3.6(-3) 


1.3(-3) 


3.1(-4) 


7.3(-5) 


14 


1.3(-3) 


4.9(-4) 


1.3(-4) 


2.4(-5) 


8.7(-6) 


8.2(-6) 


16 


1.2(-4) 


3.0(-5) 


4.7(-6) 


2.8(-6) 


2.7(-6) 


2.6(-6) 


18 


l.l(-5) 


5.8(-6) 


8.0(-6) 


8.1(-6) 


8.0(-6) 


8.0(-6) 








Outflow C 


ase 






L 


7 = 1.2 


7 = 1.4 


7 = 1.666 


7 = 2.0 


7 = 2.5 


7 = 3.0 


8 


5.4(-3) 


5.4(-3) 


5.4(-3) 


5.4(-3) 


5.4(-3) 


5.4(-3) 


10 


9.2(-4) 


9.1(-4) 


9.1(-4) 


9.1(-4) 


9.1(-4) 


9.1(-4) 


12 


1.5(-4) 


1.5(-4) 


1.5(-4) 


1.5(-4) 


1.5(-4) 


1.5(-4) 


14 


2.5(-5) 


2.7(-5) 


2.0(-5) 


2.0(-5) 


2.0(-5) 


2.0(-5) 


16 


2.3(-6) 


2.6(-6) 


2.6(-6) 


2.5(-6) 


2.5(-6) 


2.5(-6) 


18 


6.6(-6) 


3.6(-6) 


8.7(-6) 


8.7(-6) 


8.7(-6) 


8.7(-6) 



Table 1. Relative errors in -D(A) for the inflow and outflow 
cases are computed by taking the maximum relative error for 
60 contour points evaluated along the semicircle <j>. Samples 
were taken for varying L and 7, leaving v + fixed at u + = 
10" 4 and v = 0.6. We used L = 8, 10, 12, 14, 16, 18, 20 and 
7 = 1.2, 1.4, 1.666, 2.0. Relative errors were computed using 
the next value of L as the baseline. 



appear to be very similar to the shock cases with one notable difference. 
These images have a small dimple near A = to account for the eigenvalue 
there, whereas those in the shock case [H [12] were computed in integrated 
coordinates and thus have no root at the origin. 

6.4. Numerical convergence study. As in [3], we carry out a numerical 
convergence study to show that our results are accurate. We varied the 
absolute and relative error tolerances, as well as the length of the numerical 
domain [—L,L]. In Tables 1-2, we demonstrate that our choices of L = 18, 
AbsTol=le-6 and RelTol=le-8 provide accurate results. 



Appendix A. Proof of preliminary estimate: inflow case 

Our starting point is Remark 12.41 in which we observed that the first- 
order eigensystem ()34|) in variable W = (w, u — v, v) T may be converted 
by the rescaling W —> W := (w,u — v,Xv) T to a system identical to that 
of the integrated equations in the shock case; see [22J. Artificially defining 
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Inflow Case 

Abs/Rel 7 = 1.2 7 = 1.4 7 = 1.666 7 = 2.0 7 = 2.5 7 = 3.0 
10- 3 /10- & 5.4(-4) iTp) 5To(-4) 3^4) 8.6(-4) 

10- 4 /10" 6 3.1(-5) 4.6(-5) 3.4(-5) 3.3(-5) 3.3(-5) 3.2(-5) 
10- 5 /!0~ 7 2.9(-6) 3.6(-6) 3.9(-6) 6.8(-6) 2.7(-6) 2.5(-6) 
10- 6 /lQ- 8 4.6(-7) 9.9(-7) l.l(-6) 6.0(-7) 2.9(-7) 3.2(-7) 

Outflow Case 

Abs/Rel 7 = 1.2 7 = 1.4 7 = 1.666 7 = 2.0 7 = 2.5 7 = 3.0 
10 _3 /10 _J> 9.2(-4) 9.2(-4) 9.1 (-4) 9.1 (-4) 9.1 (-4) 9.2(-4) 
10- 4 /10~ 6 5.3(-5) 4.9(-5) 5.3(-5) 5.3(-5) 5.3(-5) 5.3(-5) 
10- 5 /10 -7 6.7(-5) 6.7(-5) 6.7(-5) 6.7(-5) 6.7(-5) 6.7(-5) 
10- 6 /lQ- 8 2.9(-6) 2.9(-6) 2.9(-6) 2.9(-6) 2.9(-6) 2.9(-6) 
Table 2. Relative errors in D(\) for the inflow and outflow 
cases are computed by taking the maximum relative error 
for 60 contour points evaluated along the semicircle <fi. Sam- 
ples were taken for varying the absolute and relative error 
tolerances and 7 in the ODE solver, leaving L = 18 and 
7 = 1.666, v + = 10~ 4 , and vq = 0.6 fixed. Relative errors 
were computed using the next run as the baseline. 



(u,v,v') T := W, we obtain a system 
(94a) Xv + v' -u' = 0, 

(94b) Xu + u - ^tt v = — ■ 

v" (+i V 

identical to that in the integrated shock case [3], but with boundary condi- 
tions 

(95) €>(0) = v(0) = u'(0) = 

imposed at x = 0. This new eigenvalue problem differs spectrally from ([22|) 
only at A = 0, hence spectral stability of (I22p is implied by spectral stability 
of ([M|) . Hereafter, we drop the tildes, and refer simply to u, v. 

With these coordinates, we may establish (|2.3j) by exactly the same ar- 
gument used in the shock case in [3j[T2], for completeness reproduced here. 

Lemma A.l. The following identity holds for ?R.eX > 0: 
(Be(A) + |3m(A)|) I v\u\ 2 + f \u'\ 2 

(96) < V2 [ ^-\v'\\u\ +V2 f v\u'\\u\. 

Jr+ v"i J R + 

Proof. We multiply (|94b[) by vu and integrate along x. This yields 



A / u\u\ 2 + I vu'u+ I \u'\ 2 = I -p-v'u. 



R+ 



■if 
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We get (|96l) by taking the real and imaginary parts and adding them to- 
gether, and noting that |3ie(z)| + |9fm(«)| < v2|<z|- □ 



Lemma A. 2. The following identity holds for 3?eA > 0: 
(97) 



.'|2 



25Re(A) 5 



\v\ 2 + $e{\) 



\v'\ 2 1 



+ 2 



h(v) a 7 



,./|2 



Proof. We multiply (|94b[) by t/ and integrate along x. This yields 



A 



uv + 



Mg) i ,,2 



R+ « 



~(At/ WW 

R+ « 



Using (|94ap on the right-hand side, integrating by parts, and taking the real 
part gives 



A 



uv + 



h(v) 

07+ 1 



2v 2 



w'r + Ke(A) 



,'|2 



R+ « 



The right hand side can be rewritten as 
(98) 

1 



5Re 



A 



uv + 



u'v' 



h(v) 07 



|w'| 2 + ^e(A) / 



,./|2 



R+ V 



Now we manipulate the left-hand side. Note that 



A 



uv + 



u'v' = (A + A) 



uv 



u(Xv' + v") 



u'v 



uu 



-2»e(A) 
-2SRe(A) / (Au + t/)u + 



,'|2 



Hence, by taking the real part we get 



A 



uv + 



1 -1 

U V 



.'|2 



2Ke(A) 2 



This combines with (|98p to give (|97p . 
Lemma A. 3 ([3]). For /i(v) as in (|2ip . we /lave 

/»(«) 



□ 



(99) 



sup 



1~, 7 < 7, 

1 — tu_ 



where v is the profile solution to (|18|) . 
Proof. Defining 

(100) g(v) := h(v)v~^ = -u + 0(7 - 1)6~ 7 + (a + 1), 

we have g'(v) = — 1 — 07(7 — l)-t) _7_1 < for < < v < V- = 1, hence 
the maximum of 5 on £ G [i> + ,?;_] is achieved at t) = u+. Substituting (fl~9]) 
into (jlOOp and simplifying yields (|99p . □ 
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Proof of Proposition \2.3[ Using Young's inequality twice on right-hand side 
of ([Mil together with (fMl) . we get 



(»e(A) + |9m(A)|) / v\u\ 2 + 



/|2 



<V2[ ^l\v'\\u\ + V2 [ 

JR+ V 1 JR- 



u 



V u u 



< 



< 



h(v) , /,2 + (v / 2) 2 



R+ ^ 



7+1 



AO 



h{v) l2 
-^k* 2 + e 

#r 1 1 



Assuming that < e < 1 and 



kl 2 + 



7 + l 
20 2e 



v\u'\ 2 + — 
4e 



v\u\ 



[1 — e)/2, this simplifies to 



(»e(A) + |9fm(A)|) / 6|< + (1 - e) 
1 - e 



/|2 



< 



h Wy\ 2 i 



7 + l 
20 2e 



DM 



Applying ([97|) yields 

(»e(A) + |9fm(A)|) / fi|u| 2 < 
or equivalently, 

(»e(A) + |3m(A)|) < 
Setting e = 1/(2^/7 + 1) gives (J3U). 



1 - e 2e 



(2 7 -l)e + l 



2e(l 



□ 



Appendix B. Proof of preliminary estimate: outflow case 



Similarly as in the inflow case, we can convert the eigenvalue equations 
into the integrated equations as in the shock case; see [22J. Artificially 
defining (u, v, v') T := W, we obtain a system 



(101a) 
(101b) 



Xv + v — u 



0. 



Xu + u' - ttMV - 



identical to that in the integrated shock case [3], but with boundary condi- 
tions 

(102) 



v'(0) 



X 



a — 1 



v(0), u'(0) = au'(0) 



imposed at x = 0. We shall write wq for w(0), for any function w. This new 
eigenvalue problem differs spectrally from (j22|) only at A = 0, hence spectral 
stability of (f22|) is implied by spectral stability of (jlOip . Hereafter, we drop 
the tildes, and refer simply to u, v. 
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Lemma B.l. The following identity holds for 5ieA > 0: 



(Ke(A) + |3fm(A)|) / v\u\ 2 



v x \u\ + 



i / 12 •*> i i2 

\u I + -fol^ol 



(103) 



< V2 



h(v) . ... i , I, / 1, 

\v \\u\ + I v\u \\u\ + v2|a||w ||no|. 



Proof. We multiply (|101bp by -On and integrate along x. This yields 



A 



v\u\ + 



U 



i\2 



h(v) ,_ , _ 

V U + U n Un. 



We get ()103p by taking the real and imaginary parts and adding them to- 
gether, and noting that |5Re(z)| + |9m(z)| < \pl\z\. □ 

Lemma B.2. The following inequality holds for 3fteA > 0: 



h{v) 07 
^TT + ^7+1 



v'\ 2 + Ke{X) 



/|2 



1 2 



+ 



E- V 



Av 



+ 25fteA 



(104) 



< 



1 / 1 9 ~ i i9 
|li I + Vq\Uq\ . 



Proof. We multiply (jlOlbp by v' and integrate along x. This yields 



A 



uv + 



- £7+1 



,'|2 



R- V 



R- « 



-(Av' + i/V- 



Using (llOlap on the right-hand side, integrating by parts, and taking the 
real part gives 



A / uv' + 



1 -1 

u v 



h{v) Vx_ 

g7+l + 2 ?)2 



v'| 2 +^e(A) 



/|2 



- + ■ 



R- « 



|2 



2*o 



The right hand side can be rewritten as 



(105) 



uv + 

m 

&Y+ 1 



+ 



07 



|w'| 2 + Ke(A) / 
Jm- 



,'|2 



„/ 12 



+ 



2v 



Now we manipulate the left-hand side. Note that 



uv + 



u'v' = (A + A) 



{u'v' — Xuv') 



-2Ue(X) / u'v + 2ReXu v + / u'(v' + Aw) - A-u v 
Jr- Jr~ 

-23?e(A) / (At; + t>')v + / |i/| 2 + 25fteAu vo - Au % 



Hence, by taking the real part and noting that 

!Re(2;ReAuovo — Xuqvq) = KeAKe(novo) — $smX$sm(uovo) = $te(XuoVo) 
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we get 



A / uv' + 



i -i 

U V 



u'\ 2 -2^e\ 2 \ M 2 -5fteA|uo| 2 +5fte(Au uo). 



This combines with (|105j) to give 

\v'\ 2 + ^e(X) 



h{v) 07 



\v'\ 2 \v'\ 2 



+ $te\\v \' 



2v 

u'\ 2 + $le(\uoVo). 



We get (|104|) by observing that (|102p and Young's inequality yield 



,/ 1 2 



|3Re(Aw ^ )l < \a-l\\v' v \ < \v' v \ < + v \u \ 2 . 



Here we used la — II 



I A — r>; 



tt < 1. Note that KeA > and % < 0. 



□ 



Proof of Proposition \2.3[ Using Young's inequality twice on right-hand side 
of (|103p together with (|99p. and denoting the boundary term on the right 
by h, we get 



(»e(A) + |9fm(A)|) / v\u 



< y/2 



v x \u\ 



\ u '\ 2 + t^oKP 



h(v) . ... . , .,/,,, 

|f 1i + / V\U \\U\ + ifo 



< 



< 



R- « 



7+1 



20 



/|2 



v\u\ + e / urn + 



4e 



v\u\ + A 



^§|,'| 2 + e 
£)7+l 1 1 



kl 2 + 



Here we treat the boundary term by 



JL + I' 

20 4e 



VI2 1 



i& < v^|a|K|K| < + -|a| 2 t)o|no| 2 . 

z f " 

Therefore using (jl04l) . we simply obtain from the above estimates 

(5Re(A) + |9fm(A)|) / v\u\ 2 + (1 - e) f \u'\ 2 + ]-v \uo\ 2 
Jr- Jr- 2 



< 



< 20 



h(v) . , 



J |2 



I '|2 , 



2 v 

7 r 

20 4e 



+ 



7 1 
20 4e 



■0|li| 2 + — I Q;| 2 t>o I no I 



where J\, := (^|a| 2 + 26)vq\uq\ 2 . Assuming that e + 20 < 1, this simplifies to 



(Se(A) + |3m(A)|) / £|u| 2 + -^oNI 2 < 



7 + I 
20 4e 
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Note that \a\ < < Therefore for |A| > we get |a| < 9 and 
Jb — 39vo\uq\ 2 . For sake of simplicity, choose 6 = 1/6 and e = 2/3. This 
shows that Jj, can be absorbed into the left by the term iw |^o| 2 and thus 
we get 



(5fte(A) + |9fm(A)|) / v\u\ 2 < 



7 + l 
26 4e 



v u 



3 7 + 



V u 



provided that |A| > 1/(4(9) = 3/2. 
This shows 



(»e(A) + |3m(A)|) < max{— , 3 7 + -}. 

2 8 



□ 



o 



Appendix C. Nonvanishing of D 

Working in (v,u) variables as in (|94p . the limiting eigenvalue system and 
boundary conditions take the form 



(106a) 
(106b) 



Xv + v' — u = 0, 

Xu + u v = — . 

V V 



corresponding to a pressureless gas, 7 = 0, with 

(107) (u, u, v, v )(0) = (d, 0, 0, 0), (u, u', v, v')(+oo) = (c, 0, 0, 0). 
Hereafter, we drop the tildes. 



Proof of Proposition \S '.61 Multiplying (|106b|) by vu/(l — v) and integrating 
on [0, b] C M + , we obtain 



A 



V I I 2 7 

\u\ dx + 



v 



u'udx — / v'udx - 
o l-« Jo 7o 1-w 



/o 1 - « 

Integrating the third and fourth terms by parts yields 

/- 



b u"u 



dx. 



A / -\u\ dx + 

/o l-« Jo 



V 



+ 



1 



1 — V V 1 — V 



b \u'\ 2 



u'udx 



„dx+ / v(Xv + v')dx 
o 1-u Vo 



u'-u 



1 - u 



Taking the real part, we have 
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5?e(A) 

(108) = 5Re 

where 



v 



1-v' 
u 

vu + 



u\ 2 + \v\ 2 \dx+ I g(v)\u\ 2 dx + 



9{V) 



1-v 2 



v I 1 

+ 



1 — V \ 1 — V 



+ 



1-v 



1-v 



Note that 
d 



1 



dx \ 1 — v 



{1-v)' 



v(v — 1) 



(i-v) 2 (i-v) 2 (i-v)< 



b 

1-v 



-dx 



1-v 



Thus, g(v) = and the third term on the right-hand side vanishes, leaving 



»e(A) 



1 - v 



-\u\ + \v\ I dx + 



-dx 



^ke(vu) + 
^ke(vu) + 



1 - v 
1-v 



21 



(b). 



We show finally that the right-hand side goes to zero in the limit as 
b — > oo. By Proposition 14.3^ the behavior of u, v near ±oo is governed by the 
limiting constant -coefficient systems W = ^4j_(A)VK, where W = (u,v,v') T 
and A± = A°(±oo, A). In particular, solutions W asymptotic to (1,0,0) at 
x = +oo decay exponentially in (u',v,v') and are bounded in coordinate u 
as x — > +oo. Observing that 1 — v — > 1 as x - > +oo, we thus see immediately 
that the boundary contribution at b vanishes as b — > +oo. 

Thus, in the limit as b — > +oo, 

f' +oc ( V \ f' + °° \v'\ 2 

(109) Ue(X) / — -r|u| 2 + M 2 Mx+ / = 0. 

But, for KeA > 0, this implies v! = 0, or u = constant, which, by u(0) = 1, 
implies n = 1. This reduces (|106a|) to v' = Xv, yielding the explicit solution 
v = Ce Xx . By v(0) = 0, therefore, v = for JfteA > 0. Substituting into 
(|106b|) . we obtain A = 0. It follows that there are no nontrivial solutions of 
(HMD, AMD for KeA > except at A = 0. □ 

Remark C.l. The above energy estimate is essentially identical to that used 
in |12j to treat the limiting shock case. 
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Appendix D. Nonvanishing of D k ^ 

Working in (v,u) variables as in (|94p . the limiting eigenvalue system and 
boundary conditions take the form 



(110a) 
(110b) 



Xv + v — u = 0, 

Xu + u v = — . 

v v 



corresponding to a pressureless gas, 7 = 0, with 
(111) (u, u', v, v')(-oo) = (0, 0, 0, 0), 



(112) 

In particular, 
(113) 



ff(0) 



A 



a - 1 



v(0), fi'(0) = av'(0). 



fi'(0) 



Xa 



a — 1 

Hereafter, we drop the tildes. 



u(0) = v\0)v(0) = (v - l)v v{0). 



Proof of Proposition [X51 Multiplying (jllObp by vu/(l — v) and integrating 
on [a, 0] C M - , we obtain 



A 



o 



— Inl dx + 



-u'udx — I v'udx 



o 



u"u 



1 - V J a 1 - V J a .L 1 - V 

Integrating the third and fourth terms by parts yields 



-dx. 



X 



1-v 



-\u\ dx + 



+ 



v ( 1 

1 — v \ 1 — V 




u'udx 



vu + 



\uf 
l-v 
u'u 



dx + v(Xv + v')dx 



l-v 



Taking the real part, we have 



5Re(A) 



1 - v' 



u\ 2 + \v\ 2 I dx + / g(ii)\u\ 2 dx + 



\u'\ 2 
l-v 



dx 



(114) 
where 



5Re 



vu + 



u'u 



1 



l-v 2 



v I 1 

+ 



1 — V \ 1 — V 



u 
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and the third term on the right-hand side vanishes, as shown in Section IC| 
leaving 



5Re(A) 



1-v' 



\u\ + \v\ dx + 



0^/12 
1 - V 



-dx 



$te(u'u) \v\ 
1-v 2 



A boundary analysis similar to that of Section O shows that the contri- 
bution at a on the righthand side vanishes as a — * — oo; see |12] for details. 
Thus, in the limit as a — > —oo we obtain 



Ke(A) 



l-ij [ 



f° \u'\ 2 

\u\ 2 + \v | 2 I dx + I -dx 

,1-v 



$te(vu) + 



1-v 



(0) 



< 



(1 — vo)dte(vu) 



(1 - uo)|v||«| 



(0), 



(0) 



< (1 - v ) 



Ko)l s 



where the second equality follows by (|113p and the final line by Young's 
inequality. 

Next, observe the Sobolev-type bound 



Ko)| 2 < 



v 2 

|?/(x)|<ix) < 



k" 2 



1-v 



-{x)dx 



[1 — v)(x)dx, 



together with 

(1 — v)(x)dx 



o d , 



— {x)dx = I (log?) 1 )'(x)dx = log v , 



hence /^(l — v){x)dx < nz^yz f° r v o > v *i where v* < e 2 is the unique 
solution of 

(115) v, = o- 2 /^-**) 2 

Thus, for vq > v*, 



(116) 



for e : = 



Ke(A) 



r it +M )dx + e J— ^ 

oo X 1 - « / ./-oo 1 - w 



< 0, 



(l-^o) 2 



-> — - to — ; — 1 w \ — > 0. For 5fteA > 0, this implies v! = 0, or 

I- oc ( 1 ~ v )( x ) dx 

u = constant, which, by u(— oo) = 0, implies u = 0. This reduces (jllOap 
to w' = Au, yielding the explicit solution u = Ce Xx . By u(0) = 0, therefore, 
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v = for JfteA > 0. It follows that there are no nontrivial solutions of (jllOl) . 
(fTill) for KeA > except at A = 0. 

By iteration, starting with v* ~ 0, we obtain first < e~ 2 ~ 0.14 then 
v* > e 2 /^- 14 ) 2 « .067, then «, < e 2/(i-067) 2 ra 10; then ^ > e 2/(i-.io)^ w 
.085, then t>* < e 2/(i--085) ^ 091 and ^ > g 2/(i-.09i) _ _ 889, terminating 
with « .0899. □ 

Remark D.l. Our Evans function results show that the case vq small not 
treated corresponds to the shock limit for which stability is already known by 
[12j . This suggests that a more sophisticated energy estimate combining the 
above with a boundary-layer analysis from x = back to x = L + S might 
yield nonvanishing for all 1 > Vq > 0. 



Appendix E. The characteristic limit: outflow case 

We now show stability of compressive outflow boundary layers in the 
characteristic limit v+ — * 1, by essentially the same energy estimate used in 
|18j to show stability of small-amplitude shock waves. 

As in the above section on the outflow case, we obtain a system 

(117a) A5 + 5' - v! = 0, 

(117b) Xu + u'-^v' = ^. 

identical to that in the integrated shock case [3], but with boundary condi- 
tions 

(118) v(0) = — ^— 5(0), fi'(0) = av'(0). 

a — 1 

In particular, 

(119) 5/(0) = ——5(0) = v' (0)5(0). 

a — 1 

This new eigenvalue problem differs spectrally from (|22|) only at A = 0, hence 
spectral stability of (j22|) is implied by spectral stability of (|117l) . Hereafter, 
we drop the tildes, and refer simply to u, v. 



Proof of Proposition 3.7 We note that h(v) > 0. By multiplying (|117bj) by 



both the conjugate u and 5 7+1 /h{v) and integrating along x from — oo to 0, 
we have 

f° \uuv~f +1 , f° u'uv^ +1 f° , f° u"uv~< , 

-dx + \ — — - — dx — v udx = / ———dx. 



h(v) J-oo h (v) J-oo J-oo h(v 
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Integrating the last three terms by parts and appropriately using (I117ap to 
substitute for v! in the third term gives us 



A|u| 2 ^ +1 , 

4t^ dx + 



u'uv~* +l 



h{v) 



dx + v(Xv + v')dx + 



«T|«'|2 



u'udx + 



to + 



v^u'u 



h(v) 



x=0 



We take the real part and appropriately integrate by parts to get 



(120) 
5fte(A) 



£7+1 



.2 i I |2 







dx+ I g(v)\u\ 2 dx + , , 
oo J-oo h(v) 



o fi7 



u'\ 2 dx = G(0), 



where 



yr+i \ ' 



h(v) 



+ 



h(v) 



and 



G(0) 



£7+1 



+ 



h(v) \h(v) 



\u\ 2 + $te 



f ti + 



v^u'u 



If I 



evaluated at x = 0. Here, the boundary term appearing on the righthand 
side is the only difference from the corresponding estimate appearing in the 
treatment of the shock case in [181 E]- We shall show that as — > 1, the 
boundary term G(0) is nonpositive. Observe that boundary conditions yield 



vu + 



v^u'u 



x=0 



Re(«(0)«(0)) 



1 + 



W) 



x=0 



We first note, as established in [181 E], that <?({)) > on [v+,1], under 
certain conditions including the case v+ — > 1. Straightforward computation 
gives identities: 



(121) 
(122) 



jh(v) — vh'(v) = 07(7 — 1) + t) 7+1 and 
v~ t ~ 1 v x = «7 — 
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Using (|121|) and (I122D . we abbreviate a few intermediate steps below: 

g{v) = -■ 



(7 + l)«T7»(fl) - v^ +1 h'(v) d 
h(v) 2 dv 



7W 7 h(v) — v^h'(v) 



Vx_ 

' 2 



u? ((7 + l)h(v) - vh'{v)) d 



h(v~y 



+ 



dv 



h{v) 2 
7/1(7)) — vh'(v) 



(07 - h(v)) 



aVrV 



7-1 



2h(v) 



[ 7 2 (7 + 1)^ 7+2 - 2(a + 1) 7 (7 2 " 1)^ 7+1 + (a + 1) 2 7 2 (7 - 1)* 
+ a 7 ( 7 + 2)( 7 2 - l)v - a(a + 1) 7 2 (7 2 - 1)] 



(123) 



av x v' y 1 
^3" 



[( 7 + 1)07+2 + fi7 ( 7 _ 1) (( 7 + i)o _ ( a + i) 7 )< 



2/i(t) 

+ a 7 ( 7 2 - 1)( 7 + 2)t) - 0(0 + 1)7 2 (7 2 - 1)] 



> 



av x v 



7-1 



(124) 



> 



2/i(-t)) 3 

7 2 a 3 t) :E (7 + 1) 
2/i(£) 3 u+ 



[(7 + 1)£ 7+2 + a 7 ( 7 2 - 1)( 7 + 2)0 - a(a + 1)7 2 (7 2 - 1)] 



„t +1 V 



.,7+1 \ 



a 7 / 



+ 2(7-1) M- -(7-I) 



a 7 / 



This verifies <7(£>) > as 0+ — > 1. 
Second, examine 



G(0) 



VY+1 



+ 



u 



+ 



1 + 



v^v' 

W) 



$te(v(0)u(0)) 



Applying Young's inequality to the middle term, we easily get 



G(0)< 



1 



yr+i 



+ 



v 1 



1 + 



h{v) \h(v)J V h{v) 

Now observe that I can be written as 

"7V 7 - 1 lip v^v' 



|u(0)| s 



1 



l\u(0)\ 2 . 



^7+1 
/ = TT77 - 1 



h(v) 



h{v) h(v) h 2 {v) 



' ~ h 2 (v) ' 



Using (fl"2ll and (fl~22l) . we get 



07+1 

W) 



(7 - l)-0°' — + vh'(v) 



and thus 



_ (7 - i)-pT-y + vh'(v) 

~ h(v) + 



2 



h(v) h(v) h 2 (v) 



7 ~ h 2 (v) ' 



STABILITY OF BOUNDARY LAYERS 47 

Now since h'(v) = — (7 + l)v J v' + (a + l)'yv J ~ 1 v' , as v+ — ► 1, I ~ — > 0. 
Therefore, as v+ is close to 1, G(0) < 3*'(0)|u(0)| 2 < 0. This, ^(u) > 0, and 
(1120ft give, as v+ is close enough to 1, 



/o r -57+1 I /■ 

— — \u\ 2 + \v\ 2 dx+ / 
-ooVK v ) J 7- 



r0 ^ 7 

\uTdx < 0, 



which evidently gives stability as claimed. □ 

Appendix F. Nonvanishing of D in : expansive inflow case 

For completeness, we recall the argument of [19J in the expansive inflow 
case. 

Profile equation. Note that, in the expansive inflow case, we assume 
vq < v + . Therefore we can still follow the scaling (|12p to get 

< Vq < v+ = 1. 

Then the stationary boundary layer (y, u) satisfies f)15|) with vq < v + = 1. 
Now by integrating (|16p from x to +00 with noting that £>(+oo) = 1 and 
v'(+oo) = 0, we get the profile equation 

v' = v(v - 1 + a({T 7 - 1)). 

Note that v' > 0. We now follow the same method for compressive inflow 
case to get the following eigenvalue system 

(126a) \v + v' -u' = 0, 

(126b) An + v! - (/«)' 



V 



with boundary conditions 

(127) u(0) = u(0) = 0, 

where f(v) = ^Q. 

Proof of Proposition \3.8[ Multiply the equation (|126bp by u and integrate 
along x. By integration by parts, we get 

,'|2 



A 



f°° f°° , \u\ 

/ |u| 2 <ix+ / u'u + fvu-\ — dx = 0. 

Jo Jo v 



Using (|126ap and taking the real part of the above yield 

(128) KeA / \u\ 2 + f\v\ 2 dx - - / f'\v\ 2 dx+ ] -^dx = 0. 
Jo 2 J J v 

Note that 

q(7 2 -in zv> 

which together with (|128p gives KeA < 0, the proposition is proved. □ 
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